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ABSTRACT 

We have measured the cosmic momentum power spectrum from the peculiar velocities of galaxies in 
the SFI sample. The SFI catalog contains field spiral galaxies with radial peculiar velocities derived 
from the /-band Tully-Fisher relation. As a natural measure of the large-scale peculiar velocity field, 
we use the cosmic momentum field that is defined as the peculiar velocity field weighted by local 
number of galaxies. We have shown that the momentum power spectrum can be derived from the 
density power spectrum for the constant linear biasing of galaxy formation, which makes it possible 
to estimate /3s = f2°„ 6 /&s parameter precisely where O m is the matter density parameter and &s is 
the bias factor for optical spiral galaxies. At each wavenumber k we estimate /3s (k) as the ratio of 
the measured to the derived momentum power over a wide range of scales (0.026 h~ 1 Mpc < k < 
0.157 h~ 1 Mpc) that spans the linear to the quasi-linear regimes. The estimated /3s(fc)'s have stable 
values around 0.5, which demonstrates the constancy of /3s parameter at scales down to 40 h~ 1 Mpc. 
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and the amplitude of mass fluctuation as 
The 68% confidence limits include the cosmic variance. We have also estimated 
the mass density power spectrum. For example, at k = 0.1047 /iMpc -1 (A = 60 h~ 1 Mpc) we measure 
fl} n 2 Pg(k) = (2.51^094) x 10 3 (h~ 1 Mpc) 3 , which is lower compared to the high-amplitude power 
spectra found from the previous maximum likelihood analyses of peculiar velocity samples like Mark 
III, SFI, and ENEAR. 

Subject headings: cosmology: theory — cosmology: large-scale structure of universe 



1. INTRODUCTION 

Peculiar motions of galaxies in the nearby Universe are 
very powerful tool for examining the underlying mass 
fluctuations on large scales. The galaxy peculiar veloc- 
ity directly probes the large-scale matter density field, 
and thus gives biasing information of galaxy distribu- 
tion. The advantage in using the peculiar velocity data is 
that we can explore the large-scale structures in the real 
space rather than in the redshift space, without redshift 
space distortion. On the other hand, the galaxy pecu- 
liar velocities or the absolute distances have large errors 
proportional to distance because the absolute distance 
is usually inaccurately measured compared with the red- 
shift. For this reason peculiar velocity samples are in 
general much smaller than redshift survey samples. In 
spite of this disadvantage, the peculiar velocity sample 
is very useful since the observed velocity field contains 
information on larger scales compared with the density 
field in a given survey volume. 

In the study of peculiar velocity field, homogeneous 
peculiar velocity sample with well-defined selection cri- 
teria is essential because the sparse and inhomogeneous 
sample significantly induces biases in estimating phys- 
ical parameters. Over a decade or so there has been 
major progress in the peculiar velocity studies. Ob- 
servationally, there has been a dramatic improvement 
with the completion of large Tully-Fisher (TF) and D n -a 
redshift-distance surv eys in both hemisphe res. Those are 
the Mark III catalog ijWillick et al.lll997aj, that contains 
about 3,300 spiral and elliptical galaxies, the SFI catalog 
l)Giovanelli et alJl!994[) with about 1300 late-type spiral 

Electronic address: parkc(akias.re.kr;cbp@kias.re. la- 



galaxies with / -band TF distance estim ates, and the EN- 
EAR catalog ijda Costa et al.l l2000bl) containing about 
1400 early-type galaxies with D n -a distances. Recently, 
new distance indicators such as the Type la su pernovae 
l|Riess et all 119971 iRadhurn- Smith et alJ 12004ft and the 
surface brightness flu ctuation of the early-type galaxies 
( Blak eslee et al.lll999j) have been used in the peculiar ve- 
locity study. 

Theoretically, there have been many works on the pe- 
culiar velocity fields ijStrauss fe Willickl Il995t iZaroubil 
120021 for a review). The galaxy peculiar velocities have 
been used to reconstruct the three-dimensional velocity 
field or the re al-space matter densit y field of th e local 
Universe (e.g.. Ida Costa et al.lll996l iDekel et al.lll999|) . 
and to estimate the amplitude of mass fluctuation and/or 
/3 = fi°„ 6 /6 parameter, where b is the linear bias factor. 
Among the many analysis methods that have been de- 
veloped and applied to the observational data, one of 
the major development has been the POT ENT method 
( Bertschingcr ct al. 1990t lDekel et al.ll999l) which recon- 
structs the smoothed three-dimensional peculiar veloc- 
ity field from the observed radial velocity data. As a 
new m ethod of estimating /3 parameter accurately, iParkl 
(2000, hereafter POO) has introduced the cosmic momen- 
tum field and developed the momentum power spectrum 
analysis method which does not require a smoothing of 
galaxy peculiar velocities. Other studies have applied 
the velocity correlation func tion statistic ijGorski et alJ 
Il989t iBorganiet aTll2000albi). t he maximum likelihood 
(ML) method <|Freudling et alJ 119991 ISilbcrma net alJ 
120011). the orthogonal mode expansion (OME ) method 
ijNusser fc David Il995t Ida Costa et all 119981), Wiene r 
filtering method l)Zaroubi. Hoffman, fc DekeJ j 1999) , 
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the unbiased min imal variance (UMV) estimator 
llZaroubi et al.l 12002ft. t he optimal moment expansion 
method (Wat kins et alJ |2002). and the pairwise velocity 
method (F eldman et alJl2003 ). Considerable efforts have 
been made to estimate the parameter from the compar- 
ison of mass density field derived from the observed ra- 
dial velocities with the observed galax y density field from 
the redshift survey (6-6 comparison: ISiead et al-l figOS'). 
or from the direct comparison of the observed peculiar 
velocity field with that derived from the galaxy red shift 
survey (v-v comparison; e.g., see lZaroubi et al.ll2002L ref- 
erences therein). 

In this paper, we estimate the (3 parameter and the am- 
plitude of mass fluctuation from the peculiar velocities of 
the SFI galaxy sample by measuring the momentum and 
density power spectra. The outline of this paper is as 
follows. In §2, we briefly review the cosmic momentum 
field. Description and data reduction of the SFI galaxy 
sample are given in §3. In §4, we measure the density 
and momentum power spectra from the SFI data. Mea- 
surements of (3 parameter and the amplitude of mass 
fluctuation are given in §5. We discuss our results in 
§6. Throughout this paper, in calculating the real-space 
distance from the redshift or the recession velocity, we as- 
sume a flat ACDM universe with matter and dark energy 
density parameters f2 m = 0.27 and = 0.73, respec- 
tively. For the present value of the Hubble parameter, 
we use Ho = lQOh km/s/Mpc. 

2. COSMIC MOMENTUM FIELD 

This section summarizes the cosmic momentum field. 
For complete description we refer to POO. The cosmic 
momentum field is defined as the peculiar velocity field 
v weighted by local density p/p: 



p = —v = (1 + 6)v. 
P 



(1) 



Here 6 = (p — p)/p is the dimensionless density contrast. 
The cosmic momentum p defined here has the dimension 
of velocity and is equal to v in the linear regime (\6\ <C 1). 

To compare the observed peculiar velocities of galaxies 
with cosmological models we measure correlation func- 
tion (CF) or power spectrum (PS) of the momentum 
field. Suppose the matter over-density 6(x) and the pe- 
culiar velocity v(x) fields are homogeneous and isotropic 
random fields. Given a momentum correlation tensor 
S,fj(r) — (pi(x)pj(x + r)), we define a power spectral 
tensor P[Jk) as 



d 3 r&(ry h 



(2) 



With this definition, Pfj(k) has an unit of [velocity] 2 x 
[volume], conventionally km 2 s~ 2 (ft, -1 Mpc) 3 . We call 
its trace P p (k) = i*(fe) = f d 3 x£ p (r)e lk x as the 
PS of the momentum field, where £ p (r) = £fi( r ) = 
(p(x) ■ p(x + r)) is the dot-product CF of the momen- 
tum field. 

What we directly measure in a galaxy peculiar velocity 
survey is the radial component of peculiar velocities at 
the locations of galaxies. The radial component of mo- 
mentum, the physical observable in our analysis, is the 
galaxy number-weighted quantity p r (x) = u(x)n(x)/n, 



where the radial peculiar velocity u(x) is caused by the 
total matter field and n(x)/n represents the distribution 
of galaxies. With the isotropy of the momentum field 
and the far-field approximation, POO obtained 



P p (k)^3P Pr (k). 



(3) 



The three-dimensional momentum PS can be measured 
from the PS of radial component of momentum P Pr (k). 
This property can be directly applied to the all-sky pe- 
culiar velocity survey data. By comparing the PS of the 
total momentum with those of the radial component of 
the momentum vector observed at a corner or at the cen- 
ter of the simulation cube, POO has demonstrated that 
equation (3) is actually very accurate over wide scales 
when the cosmic variance and observational uncertain- 
ties in the PS are taken into account. However, at scales 
corresponding to the fundamental modes within the sur- 
vey volume, equation (3) holds only approximately. 

If the velocity field is irrotational (curl- free), Fourier 
mode of velocity field in the linear regime becomes 



v(k) = -i(DHf)^~6(k), 



(4) 



where D(t) is the linear growth factor as in 6(k;t) — 
D(t)6(k), H = a/a is the Hubble parameter, a(t) is the 
expansion fact or, and f(fl m ,flA) = d\a.D/d\aa ~ f2^ 6 
l)Peebleslll980D . We use the Fourier transform (FT) con- 
vention defined as v(k) — (1/V) J v d 3 xv(x)e lk ' x on a 
large volume V over which they are considered to be 
periodic. The inverse FT is thus defined as v(x) — 
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V/(2tt) 3 Jd 3 kv(k)e- lli x = £ fc u(fc)e 

From equation (4), the PS of velocity field (P v ) in the 
linear regime is related with the density PS (Ps) as 



P v (k) = V(\v(k)\ 2 ) = (DHf) 2 P s (k)/k 2 



(5) 



From the FT of equation (1), p(k) — v(k) + 
J2k' 3(k')v(k — k'), the approximate expression for the 
momentum PS is given by the sum of PS of v and 6v 
fields (POO): 



P p (k) wP„(fc)4 
+ l(D 2 Hf) 2 



Psv(k) = 
d 3 k' 



(DHf) 2 P & (k)/k 
k 2 
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2 P s (k')P 5 (\k-k'\ 



(6) 



It should be noted that equation (6) is not a complete 
expression for the cosmic momentum PS because only 
the 6v term is considered with the other non-linear terms 
excluded. Furthermore, equation (6) assumes that 6(x) 
and v(x) are Gaussian with a property that the ensemble 
of odd-product like (6(x)v(x) ■ v(x + r)) is zero. The 
correct expression for the mom entum PS requires higher 
order perturbation theory fe.g.. lBernardeau et alJ l2002). 

The accuracy of equation (6) can be tested using the 
A^-body simulation data. Figure 1 shows the momentum 
and density PS directly measured from a iV-body simu- 
lation of the ACDM cosmologi cal model (open a nd filled 
circles). A particle-mesh code l|Parldll990l Il997j) is used 
to gravitationally evolve 256 3 CDM particles from z = 23 
to on a 512 3 mesh whose physical size is 614.4 /i _1 Mpc. 
The cosmological parameters used in the simulation are 
O m = 0.27, n b = 0.0463, n A = 0.73, and h = 0.72, 
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Fig. 1. — Momentum and density PS measured from the TV- 
body simulation of the ACDM universe (open and filled circles), 
in units of km 2 s _2 (/i -1 Mpc) 3 and (/i _1 Mpc) 3 , respectively. The 
bottom solid curve is the linear matter density PS (Ps) for ACDM 
cosmology as used in the Af-body simulation. The upper solid curve 
denotes the momentum PS (P p ) that is the sum of PS of v and 
8v fields (P v and Pg v ; dashed curves). The P v and Ps v have been 
derived from the linear density PS (eq. [6]) The limits of the linear 
and the quasi-linear regimes, denoted as and feQi, respectively, 
are indicated as arrows (see §5). 



and th e fitting formulae for PS given bv lEisenstein fc Hul 
( 1999) are used. The simulation is normalized so that 
er 8 = 1/6 = 0.9 at z = 0, where as is the rms mass 
fluctuation within 8 /t~ 1 M pc sphere. We use the cloud- 
in-cell (CIC) scheme (e.g., iHocknev fc Eastwood! Il98l|) 
in constructing density and momentum fields, and the 
fast Fourier transform (FFT) to estimate the PS. The 
upper solid curve which is the sum of PS of v and 5v 
fields (two dashed curves) is the momentum PS calcu- 
lated from equation (6) using the linear matter density 
PS (bottom solid curve). Comparison of the momentum 
PS measured from the TV-body simulation (open circles) 
with the upper solid curve shows that the equation (6) is 
very accurate down to scales of about 50 ft. _1 Mpc (at the 
largest scale corresponding to the simulation box size, the 
agreement is not good because of the cosmic variance). 
This demonstrates that the momentum PS can be de- 
rived from the linear density PS even in the quasi-linear 
regime. 

3. THE SFI GALAXY SAMPLE 

3.1. Description of the Sample 

The SFI galaxy catalog is an all-sky sample of pecu- 
liar velocity of galaxies. It conta ins 1289 Sbc-Sc galax- 
ies with /-band TF distances ijGiovanelli et al J 11994. 
Il997albt IHavnes et alJ 1999a bl. The galaxies in the 
catalog have inclination i > 45° and Galactic latitudes 
|6| > 10°. Its survey depth is czlg = 7500 km/s, 
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Fig. 2. — Redshift (cz^q) versus angular diameter (d m ) of 1289 
galaxies in the SFI catalog. The d m is obtained from a relation 
dm/0,23.5 = 2.78 — 1.03 log(a/6), where 023.5 is the isophotal angu- 
lar radius measured at fii = 23.5 mag/arcsec 2 level and a/b is the 
axial ratio. The region enclosed by solid lines represents the SFI 
main galaxy selection criteria, while the dashed line denotes our 
angular diameter limit applied to the second and the third redshift 
shells. 



where czlg is the galaxy recession velocity with respect 
to the Local Group (LG). This sample is a combina- 
tion of the SFI main g alaxy sample (S > —45°) and 
IMathewson et alJ ij 1992. hereafter MAT) sample of which 
the magnitudes and rotational velocities were converted 
to the SFI system. 

The peculiar velocities of spiral galaxies in the SFI 
sample have been derived from the /-band TF relation. 
We consider the TF relation between the absolute mag- 
nitude Mj and the velocity width parameter rj, 



Mj = nij — 51ogr = otf + ^tf 7 ?, 



(7) 



where rj = logw — 2.5 and w is the circular velocity line- 
width of a galaxy in unit of km/s. The absolute mag- 
nitude Mj is defined as the apparent magnitude when a 
galaxy is located at r = 1 km/s. For the slope and the 
zero-point of the TF relation, we adopt axF = —21.10 
and &tf = -7.94 of the in verse Tully-Fisher (ITF) re- 
lation that is determined bv lGiovanelli et aD <)1997a|) for 
24 clusters in the SCI sample. The reason for using the 
ITF relation is worth mentioning. The forward TF rela- 
tion is obtained by regressing the apparent magnitudes 
over the line-width. It can be biased due to the im- 
posed selection limits on magnitude, angular diameter, 
and circular velocity line-width. On the other hand, 
the ITF relation is obtained by fitting the line-width 
as a function of the apparent magnitude. It avoids the 
selection bias i f the sample selection is independent of 
the l i ne-width llStrauss fc WillicldlT995l iFreudling et alJ 
119951 iBorgani et alJl2000afl . We perform our analysis us- 
ing the peculiar veloci ties inferred from the ITF relation. 
IBorgani et alJ l)2000a|) analyzes the SFI peculiar velocity 
data in redshift space to avoid the possible bias that can 
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SGX (km/s) SGX (km/s) 

Fig. 3. — (a) Hammer- Aitoff sky projection in the Galactic co- 
ordinate as seen in the LG-frame, of the 661 SFI galaxies. Open 
circles indicate the infall while crosses outflow. The size of symbols 
scales with the velocity amplitude as shown in the bottom-right of 
the figure, (b) & (c) Distribution of the SFI galaxies in the su- 
pergalactic (SG) coordinates before and after MB correction, re- 
spectively. The positions and peculiar velocities are shown in the 
SGX-SGY plane with thickness of AZ < 2000 km/s. Outflowing 
galaxies are denoted as dots with solid segments that scale with 
the velocity amplitude, while infalling galaxies as open circles with 
dotted segments. The large circle in panel (c) represents a region 
enclosed by a 40 h _1 Mpc sphere. 



enter into the measurement of distances or peculiar ve- 
locities of galaxies. However, we make our analysis in 
real space to avoid the effects from the redshift space 
distortion. 

The SFI galaxy sample, by desire, is angular-diameter 
limited. The sample was made to obtain roughly the 
same number of galaxies in three different redshift shells. 
Diameter limits imposed on the SFI main galaxy sample 
on each redshift shell are 2 .'5 < d rn < 5' at cz LG < 3000 
km/s, I'.Q < d m < 5' at 3000 < cz LG < 5000 km/s, 
and 1.'3 < d m < 5' at 5000 < cz LG < 7500 km/s. 
Here d m is the angular diameter defined as the UGC 
blue major axis o r the analogous param eter in the ESO- 
Uppsala Survey ijGiovanelli et aljfl994|) . The distribu- 
tion of SFI galaxies on a cz^Q-dm plane is shown in Fig- 
ure 2. The galaxies actually do not strictly follow the 
selection criteria. The angular diameter of each galaxy is 
estimated from a relation d TO /a23. 5 = 2.78—1.03 log(a/6), 
where 023.5 is the /-band isophotal radius measured at 
/ij = 23.5 mag/arcsec 2 level and a lb is the axial ratio 
of a galaxy ijGiovanelli et alJl!994|) . The angular diam- 
eter limit of the MAT galaxy sample is 1.'7. For com- 
pleteness of the SFI sample, we exclude galaxies with 
diameter smaller than 1'.7 on the second and the third 
redshift shells. A small fraction of galaxies with small 
line-widths (logw < 2.25) have also been discarded be- 
cause of the large fractional errors in the measurement of 
their widths and thus the la rge uncertainties in t he pe- 
culiar velocity measurement (|da Costa et al.l l"l996). The 



total number of galaxies we use in our analysis is 661. 
Figure 3 shows the positions and peculiar velocities of 
the 661 SFI galaxies on the sky and the supergalactic 
plane. The supergalactic coordinate system is a Carte- 
sian coordinate system with the SGZ-axis pointing to the 
Galactic coordinate I = 47?37, b = 6?32 and the SGX - 
axis pointing to I = 137?29, 5 = 0° fe.g.. lPeebleJll99l| . 

3.2. Bias Correction 

Real data are often selected by magnitude or diame- 
ter limits, and/or by redshift limit, which cause a bias 
in the calibration of the TF relation. This is particu- 
larly important for the SFI sample because the adopted 
redshift-dependent selection criteria imply that near each 
redshift limit outflowing galaxies are preferentially ex- 
cluded from the sample. As shown in Figure 36 this effect 
becomes dominant at large distances, leading to a spuri- 
ous systematic infall of galaxies at th e outer edge of the 
survey volume (|da Costa et alJ ll996). In earlier studies 
with the SFI data, the biases wer e estimated using the 
numerical Monte-Carlo technique llFreudling et aflll995t 
Ida Costa et a,l.lll996t iBranchini et alJl2001l) or using the 
semi-analytic approach ijFreudling et al.1 1999). In this 
paper, we incorporate all the selection effects into two se- 
lection functions: a radial function 4>{r) for the redshift- 
dependent selection, and an angular function ip(b) for the 
Galactic extinction that may reduce the galaxy number 
density at the low Galactic latitude. For the radial se- 
lection function we use a method that is equivalent to 
the V/VniM method for computing the luminosity func- 
tion llSchmidtlll968tl. The radial selection function 4> {r) 
is given bv l|Park et alJll99l iStrauss fc WllTiclll995l) 

'M = J- V -J-, (8) 



0, 



where Q s is the solid angle covered by the survey and 
d maJi i is the maximum distance the i-th galaxy can have 
while satisfying the selection criteria. 

Malmquist bias is caused by the random error in the 
galaxy distances estimate d by the distance estim ators 
like TF or D n -a relations <|Lvnden-Bell et al.lll988h . At 
a given distance there are more galaxies that are per- 
turbed from the far side than from the near side. The 
number of galaxies that are randomly moved to that dis- 
tance also depends on the actual distribution of galaxies. 
Biases in the inferred distance or the peculiar velocity 
occur because of the volume effect and galaxy number 
density variation along the line of sight. We call the 
combination of the two effects the Malmquist bias (here- 
after MB). To correct for MB, we have derived a density 
distribution of the local uni verse using the IRAS PSCz 
0.6 Jy flux- limited catalog ijSaunders et alJ 120001. We 
assume the linear theory of gravitational instability to 
estimate the peculiar velocity of each PSCz galaxy, and 
obtain th e real-space density fie ld using the method de- 
scribed in Branchini et al. ( 1999). In deriving the density 
field we assume (3 = 0.5. Given a TF distance (d) of each 

;alaxy, we obtain M B-corrected distance using a relation 

Strauss fc Willicldll995l 

EMd) _ Jo" rMr)Hr) exp (~[ln(r/d)] 2 /2A 2 ) dr 

r o °°r 2 n(r)0(r)exp(-[ln(r/d)] 2 /2A 2 )dr' 1 ' 

Here n(r) is the number density of galaxies along the 
line-of-sight, r denotes unbiased true distance, and A = 
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Fig. 4. — Angular (upper) and radial (bottom) selection functions 
of the 661 galaxies in the SFI sample. For the angular selection 
function we estimate the number density of galaxies at Galactic 
latitude bins with Ab = 10° . The estimated number densities 
have been fitted with a function ip(b) = i/>olO Q ' 1-csc I 6 " , where 
V>0 = 81.3 (75.2) and a = 0.113 (0.084) for the northern (south- 
ern) hemisphere. Radial selection functions before and after MB 
correction, and that further smoothed with A = cr TF lnlO/5 are 
shown as dotted, solid, and dashed curves, respectively. 



cttf In 10/5 is a measure of the fractional distance uncer- 
tainty, where cttf is the error of the TF relation in unit 
of magnitude. In this paper, we assume cttf = 0.36 
ijGiovanelli et alJ fl997a|) . The (f>{r) in equation (9) is 
the real-space radial selection function derived from the 
raw SFI sample {not corrected for MB; dotted curve in 
Fig. 4). Given the fully corrected distance estimate, 
d c = E(r\d), the radial component of the peculiar veloc- 
ity is obtained as u = cz-^Q — dc. The computed distances 
and peculiar velocities are used to measure the cosmic 
density and momentum PS (§4.3). 

Figure 4 shows the angular and radial selection func- 
tions of the 661 galaxies in the SFI catalog. The angular 
selection function (ift) is obtained by fitting the galaxy 
number density estimates at Galactic latitude bins to a 
function tp(b) — V'oe"*- 1-080 > b *K We calculate radial selec- 
tion functions (<fi) from the 661 galaxies before and after 
MB correction (dotted and solid curves, respectively), 
and MB-corrected selection function smoothed with a 
Gaussian filter of width A (dashed curve). The final se- 
lection function is given by s(r, b) = 4>(r)^p{b). 

3.3. Survey Window Function 

Although galaxy positions in redshift space are mea- 
sured very accurately, the error in the absolute distances 
and the peculiar velocities of galaxies monotonically in- 
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Fig. 5. — (Upper) Effective weight of a galaxy at distance r (solid 
curve), shown as the product of the inverse of the radial selection 
function (<j>~ 1 (r); dotted) and the survey window function (W(r); 
dashed curve) for <ro = 200 km/s and (Ttp = 0.36. Here <j>(r) is 
normalized so that (f>(0) = 1. (Bottom) PS of the survey window 
function |W(fc)| 2 . The survey region is limited to a volume within 
a sphere of r m ax = 40 fe-iMpc and |6| > 10°. 



creases radially outward, and dominates the signal be- 
yond a certain distance. The expected distance error 
of a galaxy at distance r (in /i Mpc) is approximately 
u(r) = 100r(e A — 1) km/s. To reduce the effects of the 
large distance error on the measured PS, we use a survey 
window function defined as 
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With the smaller cto, the effective survey depth in our 
analysis becomes shallower. Figure 5 shows the inverse 
of radial selection function _1 (r) and the window func- 
tion W(r) (dotted and dashed curves). The product 
0~ 1 (r)W(r), the solid curve in Figure 5a, is the effective 
weight of a galaxy at distance r (see eq. [12] in §4.1). 
We adopt cto = 200 km/s and limit the sample depth to 
r nra = 40 /i _1 Mpc to exclude distant galaxies with large 



() 



effective weights. The PS of the survey window function 
shown in Figure 5b is a monotonically and rapidly de- 
creasing function, and is negligible at k > 0.1 ZiMpc -1 . 
That is, the measured PS has correlations only within 
about Afc = 0.1 /iMpc -1 . 



4. POWER SPECTRUM ESTIMATION 
4.1. Method 

In this section we describe the method of estimating 
the PS of the momentum field P p {k) from the observed 
radial peculiar velocity data. We use the PS estimatio n 
method developed by POO (see also iPark et alJ [1994). 
The method uses a direct FT to calculate the Fourier 
modes of density and momentum fields. In this paper we 
slightly modify the direct FT method so that it can be 
applied to a case for an arbitrary window function. 

Suppose we have N galaxies with positions Xj and ra- 
dial peculiar velocities u(xj), with j = 1, . . . ,N. Owing 
to the variation of selection, each galaxy has a different 
statistical weight that is given by the inverse of the selec- 
tion function, wj — s^ 1 (xj), where the sample selection 
function s(xj) is defined to be the fraction of galaxies at 
position Xj expected to be observable by the survey. If 
the galaxies represent mass, the density contrast is re- 
lated to the galaxy number density n(x) as in 



1 + S(x) 



where 5^ 3 ' (x — Xj) is the Dirac delta function, and V is 
the survey volume. 

It is quite reasonable to assume that (v) = for a suf- 
ficiently large volume. We estimate the mean peculiar 
velocity v using a bulk flow model for the radial pecu- 
liar velocities within r max . The mean velocity Vq is ob- 
tained by minimizing \ 2 — J2j ^j( u j ~~ v o ■ ^j) 2 > where 
Wj = W(xj) and Uj = u(xj), and Xj is an unit vec - 
tor pointing to the j-th galaxy l)da Costa et alJfeOQOal) . 
Then, we obtain a new radial peculiar velocity of each 
galaxy by subtracting a line-of-sight component of the 
mean peculiar velocity from the galaxy's observed radial 
peculiar velocity: v r {xj) — u(xj) — Vq ■ Xj. 

The FT of the observed radial momentum field is 



p r (k) = 1 J d 3 xW(x) [(1 + S{x)) v(x) ■ x] e ik x 
1 N 



(12) 



Likewise, the FT of S(x) is given by S(k) = 



ik X 



W(k), where W{k) is the 



FT of the survey window function. The quantities with 
a caret, p r {k) and <5(fc), are the momentum and density 
Fourier modes that are convolved with the survey win- 
dow function. The ensemble average of \p r (k)\ 2 is given 



by 

(XjWjf (\ Pr (k)\ 2 ) 

= ( WjWjVrje 



k X, \ ^ yy - 



ikx„ 



j(cells) m(cclls) 



(13) 

where j). In the second equality, we divide the 

survey volume V into infinitesimal cells with occupation 
number rij = or 1. The ensemble averaged quantity 
(njn m v r jV rm ) is related with the radial momentum CF 



{jljfljYiVrj VrtTi) 



TlSjd 3 Xj^p r (0) 
-2 

ft SjSjYi 



d Xjd x rr i^ i p r {\ [ x j 



if j = m (14) 
if j ^ to ' 



where Sj — s(xj), and we use a property (rij) — 
nsjd 3 Xj. Similarly, for density PS we need to 
consider (rijn m ). The right-hand side of equation 
(14) becomes ^ . iy|W?^ Pr . (0) for j = to case, and 

(E^) 2 £fc'|W(fc')| 2 (M fc - fc ')| 2 ) for j + m case. 
The final formulae for the density and momentum field 
PS are 



P S (k) w V 



and 



P p (fc) w 3^ 



\5(k)f 



(\Mk)\ 2 )- 



(E 7 ^) 2 



k 



E,w 2 w| 



(15) 

Ei w ( fc )i 2 ' 

k 

(16) 



where the radial momentum CF at the zero lag can be 
estimated as 



(17) 



In equations (15) and (16), we assume that PS we are 
measuring are slowly varying functions, and approximate 
that the window-convolved PS (P) can be separated into 
the true PS part (P) and the window-related part as 

P(k) Pa P(fc)Efc |W(fc)| 2 . The J2k l w ( fe )l 2 in the de- 
nominator is the overall power correction that compen- 
sates for power loss due to the finite window function 
that decreases radially outward. 

4.2. Momentum Power Spectrum from Mock Surveys 

We use the iV-body simulation data in §2 to make mock 
SFI observations. To mimic the LG environment, we 
select an observer as a particle with total velocity 600- 
700 km/s on the over-dense region with < S < 1. From 
the chosen observer we randomly select the particles with 
probability given by s(r, b). The radial selection is made 
at the true distance using the radial selection function 
derived from the MB-corrected SFI galaxy distribution 
(solid curve in Fig. 4). After the particle selection, we 
perturb the particle distances with Gaussian errors of A. 
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Fig. 6. — Momentum and density PS measured from 500 SFI 
mock observations. Shown are the median values (open and filled 
circles) and 68% confidence limits of the PS in the ACDM model. 
The solid curves are the true PS obtained from the JV-body simu- 
lation data. 



By applying redshift cut at czlg = 7500 km/s we finally 
obtain 661 angular positions and recession and peculiar 
velocities of 'galaxy' particles for each mock observation. 
For PS measurement, we use equation (9) to obtain the 
MB-corrected galaxy distances and peculiar velocities. 
Here we have smoothed the iV-body simulation data to 
obtain n(r) by applying a Gaussian filter with a = 5 
/i _1 Mpc, and used the radial selection function smoothed 
with A (dashed curve in Fig. 4). 

The PS of momentum and density fields measured from 
the five hundred SFI mock observations are shown in Fig- 
ure 6. For each mock data the number of galaxies within 
r max = 40 /i Mpc varies from 300 to 400. Shown are 
the median values (open and filled circles) and 68% con- 
fidence limits of P p {k) (i.e., 3P Pr (k)) and Pg(k) of the 
500 SFI mock PS. Note that data points in the PS are 
correlated with each other over approximately four neigh- 
boring points (Ak ~ 0.1 WVlpc -1 ). The solid curves are 
the true PS of momentum (upper) and density (bottom) 
fields of the ACDM universe measured from the ./V-body 
simulation (same as open and filled circles in Fig. 1). 
The damping of powers at large scales (or small k) occurs 
due to finiteness of the survey volume. Similar damping 
is also seen at high k, which is due to the smoothing effect 
of the peculiar velocity and distance errors. We use the 
ratio of the true PS to the PS of mock surveys as the cor- 
rection factor when we measure the PS for the observed 
SFI sample. At each fc, the raw power measured from 
the observed data is multiplied with the corresponding 
the correction factor. Then, the systematic effects in the 
raw powers are removed, and the PS with correct am- 
plitudes can be restored. Those factors also reduce any 
residual biases. This method of correction for the sys- 
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Fig. 7. — Momentum and density PS measured from the SFI 
galaxy sample (filled circles). The velocity PS measured from 
the Mark III catalog with POTENT method is shown as stars 
IKolatt &: D ckcl 1997), and the momentum PS measured by POO 
for the MAT data as diamonds. Open triangles denote the decorre- 
lated real-space galaxy density PS me asured from the SDSS galaxy 
redshift survey ITegmark et al. 2004). Open sq uares represent the 
galaxy density PS in the real space calculated by Park et al.[ 119941) 
for the 101 h~ 1 Mpc deep subsample of the CfA survey. 



tematic effects in the observed galaxy PS b y usiiu 
sur veys has been develop ed bv IPark et al J (1992, 
and lVoeelev et alJ l|1992|) . 



mock 
1994) 



4.3. Application to the SFI Data 

We analyze the SFI catalog using the same method as 
applied to the mock SFI catalogs. The number of galax- 
ies used in the analysis is 370 within r max = 40 ft, _1 Mpc. 
The final PS of the observed SFI galaxies obtained by 
multiplying the correction factor to the measured power 
spectrum at each k are shown in Figure 7 (filled cir- 
cles). The PS are measured in the true distance space 
rather than in the redshift space. The 68% uncertainty 
limit at each wavenumber is estimated from the five 
hundred mock surveys in the ACDM model. The Pois- 
son noise powers for the density and momentum PS are 
523 (fr^Mpc) 3 and 7.5 x 10 8 kmV^/i^Mpc) 3 , respec- 
tively. F or comparisons, we plo t the velocity PS mea- 
sured by IKolatt fc Dekel ((1997) who applied the PO- 
TENT method to the Mark III catalog (stars), and the 
momentum PS measured by POO for the MAT data (di- 
amonds). The velocity field is close to linear at scales 
k < 0.07 /iMpc -1 , where the momentum PS should be 
equal to the velocity P S (POO). We find that the velocity 
PS of lKolatt fc Dekel is in good agreement with the SFI 
momentum PS only at the first two wavenumbers, while 
the MAT momentum PS is consistent with the SFI PS 
at all wavenumbers. 

We also plot the real-space galaxy density PS for 



the Sloan Dig i tal S ky Survey (SDSS; open triangles; 
iTegmark et all I2004D and the Center for Astrophysics 
(CfA; open squares; iPark et aT] Il994|) galaxy samples. 
The measured SFI density PS is in good agreement with 
the SDSS and the CfA density PS even though the size 
of the SFI survey volume and the number of galaxies an- 
alyzed are much smaller. Compared to the SDSS and 
the CfA PS, the SFI density PS has relatively higher 
power at scales around 40 h~ 1 Mpc. The high power at 
those scales, considered to be a sampling effect due to 
the small survey volume, induces somewhat large value 
°f °8,s = 34, the rms mass fluctuation of the spi- 

ral galaxies within 8 /i _1 Mpc sphere. The subscript S 
denotes the optical spiral galaxies. The 68% uncertainty 
limits are from the 500 mock surveys. 

5. ESTIMATING THE f3 PARAMETERS 

Until now, we have assumed that the galaxy distribu- 
tion represents the mass field. If this is not true, the 
formulae containing the S term should be modified be- 
cause the observed momentum is now (1 + S g )v. Sup- 
pose galaxies are linearly biased with respect to mass by 
a constant factor, i.e., S g = &<5, where b is the linear bias 
factor of galaxy distribution. In the linear regime the PS 
of velocity field in equation (5) becomes 

P v (k) = (H Q (3) 2 P Sg (k)/k 2 , (18) 

where (i — f{VL m ,VL[C)/b ~ Q^/b and Ps g is the galaxy 
density PS. POO has proved that only the overall ampli- 
tude of the momentum PS given by equation (6) is scaled 
by a factor (3 2 — (f /b) 2 due to biasing if /3 is independent 
of k. 

Since we calculate both density and momentum PS 
from the same sample, it is possible to measure the /3 
parameter more accurately. We estimate /3s = ^mV^s 
parameter from the measured SFI PS. As observed by 
POO, there is a fair amount of correlation between the es- 
timated P$ and Pp. This is because density and momen- 
tum PS measured from the same sample tend to fluctuate 
statistically in a similar way. The correlation makes their 
ratios less uncertain, making the estimated /3 parameter 
more accurate. Furthermore, we note that the momen- 
tum PS can be derived from the observed density PS 
even in the quasi-linear regime (§1, Fig. 1). We estimate 
Ps from the following formula 



A [h" 1 Mpc] 
100 



/%(*) = 



pobs( 



(19) 



P p d -(fc) ' 

Here the derived momentum PS, Pp er , is calculated by 
equation (6) with the density PS replaced by Pg hs and 
D = f = 1. During the numerical integration, we inter- 
polate the density PS between the measured data points, 
and for small k < 0.02 hMpc^ 1 we extrapolate the PS 
using the ACDM linear density PS that has been scaled 
to match the observed PS. For large k > 1 hMpc^ 1 we 
also extrapolate the PS by a power-law function that fits 
the observed PS at high k. Figure 8 shows the estimated 
Ps(k) at each wavenumber k with 68% limits including 
the cosmic variance (filled and open circles with error 
bars). The uncertainties have been determined from the 
distribution of /3s parameters of the 500 mock observa- 
tions, where the PS of each mock observation are scaled 
to give the observed SFI PS on average. 
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Fig. 8. — Ps = ^m 6 /^S parameters estimated from equa- 
tion (19) (filled and open circles), together with /3l sCO = 
(fc/#o)[-Pp bs (fc)/-P,5 0bs (fc)] 1/2 (stars). The uncertainties represents 
the 68% confidence limits which include the cosmic variance. The 
first six data points (filled circles) have been used to obtain the 
weighted average of /3s parameters (see Table 1). 



We determine limits of the linear and quasi-linear 
regimes (fcj, and fcQi) by comparing the true momentum 
PS with the velocity and momentum PS derived from the 
linear density PS (see Fig. 1). The k^ (&ql) is defined 
as the scale at which the derived velocity (momentum) 
PS deviates from the true momentum PS (open circles in 
Fig. 1) by 20% in power. The two limits are k L = 0.070 
hMpc^ 1 and kq^ — 0.161 /iMpc -1 , as denoted by arrows 
in Figure 1, and correspond to wavelengths of 90 and 
39 h~ 1 Mpc, respectively. Table 1 lists the /3 parameters 
measured from the SFI galaxy sample in the range 0.0262 
/iMpc" 1 < k < 0.1571 hMpc' 1 (six filled circles in Fig. 
8). The weighted average is /3 S = 0-489±g;g|g. We call 
(3l,s as j3 parameter calculated by using equation (18) in 
the linear scales k < ki,- They are also listed in Table 
1, with a weighted average (3l,s — 0.55l ; 14 . In estimat- 
ing LjS , we include (3 at 80 /i _1 Mpc scale (k = 0.0785 
hMpc^ 1 ) where its uncertainty dominates over the 20% 
difference in power between the velocity and momentum 
PS. The Pl,s is very uncertain because the information 
contained in the survey volume is not enough to constrain 
the PS to an accurate value at large scales. It starts to 
deviate from /3s (k) significantly at k > k^ where the lin- 
ear gravitational instability theory no longer holds (see 
the stars in Fig. 8). In Figure 9 the observed SFI mo- 
mentum PS (upper open circles) is compared with that 
derived from the observed density PS (diamonds) when 
/3 = 0.489. The momentum PS prediction is accurate on 
scales down to 40 /i _1 Mpc. 

Cosmological parameters derived from the momen- 
tum PS analysis are summarized in Table 2, where /3s, 
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TABLE 1 

Ps = ^m 6 /bs and Vll^Ps Measured from the SFI Peculiar Velocity Sample 



k (hMpc" 1 ) A (ft-J-Mpc) PlM^ /%(*>) ^m 2 -P,5(fc) (h^Mpc) 3 



0.0262 


240 


42+ - 66 




4.8111; oo x 103 


0.0524 


120 


n 4Q+0-30 


u.^oa_ 147 


3.50+.^ x 10 3 


0.0785 


SO 


o w+o- 19 

U - OM -0.13 


o.4os±g:ig| 


3.nti:§| x io 3 


0.1047 


60 




502+°' 103 


2.51±°,: 9 94 x 10 3 


0.1309 


48 




469+ 121 


i.9i±J: 7 g x 103 


0.1517 


40 




n 495+0.096 


1.65t° : £> x io 3 



Average 0.55±o;?4 0.489± o ;' 



Note. - ^ iiS (fe) = (k/H )[P° ba (k)/P^ a (k)]^ 2 
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Fig. 9. — (Upper) Comparison of the observed momentum 
PS (open circles) for the SFI galaxy sample with the momen- 
tum PS (Pp CT ; diamonds) derived from the observed galaxy den- 
sity PS (Pg ha ) for /3 S = 0.489. The crosses are the velocity PS 
given by P v (k) = {H p s ) 2 P^ ha (k) /k 2 for the SFI sample. (Bot- 
tom) Derived matter density PS given by Q,]^P s (k) = PsPg bB (k) 

(circles). For compar isons, the matter density PS measured by 
IKolatt He De kcl ( 1997) for the Mark III sample (open squares) and 
by POO for the MAT sample (open triangles) are plotted. The solid 
curves are the momentum and matter density PS calculated from 
the TV-body simulation data of the ACDM universe. 



TABLE 2 

Cosmological Parameters Derived from the Momentum 
Power Spectrum Analysis 



Parameters 



Estimated Values 



P measured from SFI sample 


Ps 


- o 49+ 08 


Matter density 




- 30+ 09 


Rms mass fluctuation of spiral galaxies t 


C8,S 


- 1 14+0.33 


Amplitude of mass fluctuation 




- 56+ 27 

— u.oo_ 21 



Note. — f The rms mass fluctuation of the spiral galaxies 
within 8 /i _1 Mpc sphere obtained by integrating the observed 
SFI density PS 



tt m /b 5 s /3 (= 



cr 8i s, and ct 8 0^ 6 (= ct 8 , s /3s) are 
listed with 68% confidence limits. Here the bias fac- 
tor for optical spiral galaxies is defined to be 6s = 
cs.s/cs- Table 3 compares our 0s and osQ^ mea- 
surements with the previous results obtained from spiral 
galaxy samples (SFI, MAT, and Mark III spirals). Com- 
pared with the previous studies with the SFI da t a, our 
0S estimate is between th ose of iBranchini et all 1)200 ID 
and Ida Costa et all l)1998l) . and our gs^^; 6 between the 
large and the s mall va lues of Frcudling et al. ( 199^) and 
iBorgani et alJ l)2000a|) . By assuming that the bias factor 

1 /2 

is a function of k, given by bs(k) = [P$ g (k)/Ps(k)\ , 
we can infer the PS of mass fluctuation ^l^Ps(k) — 
0g(k)Ps g (k), by multiplying the observed density PS by 
/3|(fc) at each k (Table 1). In Figure 9 we plot the PS 
of mass fluctuation given by fl^Ps(k) — PsP$ hs (k) with 
Ps = 0.489 (filled and open circles in the bottom panel). 
For comparisons, t he matter density po wers measured 
from the Mark III (Kolatt & Dekel 1997; open squares) 
and the MAT (POO; open triangles) samples are also 
shown. The MAT mass fluctuation powers are in good 
agreement with our estimates while those of Mark III 
sample are higher than ours on 60-100 ft,~ 1 Mpc scales. 

To assess the effects of MB correction on the estima- 
tion, we have measured the momentum and density PS 
using the SFI data that are not corrected for MB. The 
measured density PS has amplitudes very similar to that 
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TABLE 3 

Comparison of Bg and crs^™ 6 Derived from the Spiral Galaxy Samples 



3a 
MS 




Data 


Method 


Ref. 


n 4Q+ 08 


5fi+ - 27 
U.OD_q 21 


SFI 


Mom. PS 


This Sti 


0.6 ±0.1 




SFI, IRAS 1.2 Jy 


OME 


1 


0.42 ± 0.04 




bl 1, 1HAS PbCz 


T 7"T7>T TV T l~\ T~\ 


2 




0.82 ±0.12 


SFI 


ML 


3 




(0.3±0.1)[r/0.2] 1 /2 


SFI 


CF 


i 




0.63 ±0.08 


SFI 


ML 


5 


51+° 13 


0.56 ±0.21 


MAT 


Mom. PS 


6 


0.6 ±0.125 




MAT 


ROBUST 


7 


0.50 ±0.10 




Mark III (S), IRAS 1.2 Jy 


OME 


8 


0.49 ±0.07 




Mark III (S), IRAS 1.2 Jy 


VELMOD 


9 


0.60« : » 




Mark III (S) 


Vrms 


10 



References . (1 ) Ida Costa et alJ il993) : pi IBranchini et alJ JgOOlT ) - (3) 

IFreudling eTai] il999l): f4llBorgani et alJ J2000afl: f5 l |Silberman et alJ J2001D: (6llPark1 
a2000n : (71IRauzv fc HendrvN2000D : (8l lDavis etafl Jl996H : f9l IWillick et all (ll997bH : 
(101 IPadilla fc Lambasl (119991) 

Note. — P is the shape parameter of the CDM models. 



from the SFI data with MB correction, but the measured 
momentum PS has lower amplitudes at all scales. The 
(3 parameter, estimated in the same way as described 
above, decreases to 0.393, which is about 2a below our 
estimation /3s — 0.489. There are systematic effects on 
(3 estimation from errors in peculiar velocities and dis- 
tances of galaxies. To measure the size of systematic 
effects, we have calculated /3 parameters using the mock 
catalogs with true distances and peculiar velocities. The 
estimated 68% uncertainty limits are (—0.13, ±0.24) and 
(-0.048, ±0.053) for f3 L<s and /3 S , respectively. They are 
very similar to but slightly smaller than those for the 
realistic SFI mock catalogs. This indicates that the un- 
certainty of j3 for the SFI sample is still dominated by 
the statistical effect rather than systematic effects. 

6. DISCUSSION 

We have measured the momentum and density PS 
from the SFI peculiar velocity data, and estimated /3s = 
£l®f/bs parameter and the amplitude of mass fluctuation 
crg^m- Our method is self-consistent because only the 
SFI peculiar velocity sample is used without relying on 
other external velocity or density fields. By noting that 
the momentum PS is accurately related to the galaxy 
density PS up to the quasi-linear regime (§2, Fig. 1) 
and that the density and momentum PS measured from 
the same sample fluctuate in a similar way, we have es- 
timated (3 parameter over a wide range of wavenumber 
space. We have determined limits of the linear and the 
quasi-linear regimes as k L = 0.070 hMpc^ 1 and Uq L = 
0.161 /iMpc^ 1 , and over the ranges k < k L and k < kQL 
we have estimated /3l,s and fls(k), respectively. Our /3s 
estimation gives stable and consistent values around 0.5 
at scales from 240 /i _1 Mpc to 40 ft. _1 Mpc (Fig. 8 and Ta- 
ble 1), with an average /3s = 0.491° 05) which translates 

into the matter density parameter f2 m = O.SO+QQg&g^ 3 . 
Our measurement is consistent with (3 parameter esti- 
mates from other studies with spiral galaxy samples, and 



is more accurate. Especially, the derived parameters are 
very similar to those from POO who has applied the mo- 
mentum PS analysis method to the MAT sample. 

Table 4 summarizes /3 or mass fluctuation measure- 
ments of other studies that use different peculiar ve- 
locity samples (ellipticals, spirals plus ellipticals, and 
Type la supernovae). In /3 estimation, most of studies 
depend on the external galaxy distribution information 
such as IRAS 1.2 Jy and PSCz catalogs for comparing 
the observed peculiar velocities with that derived from 
the galaxy redshift survey data [v-v comparison), and 
give /3 parameters around 0.5. However, the POTENT 
method, one of the S-S comparison methods, usually 
gives much higher valu e of (3 parameter ll Kolatt fc Dckcl 
ll997t|Sigad et al.ll998|1 . Compared with Kolatt fc Dekel 
(1997), our momentum PS is consistent with their ve- 
locity PS in the linear regime (Fig. 7). Despite this 
fact, their estimate of the (3 parameter is much h igher 
than ours. This is because iKolatt fc Dekell l|1997D have 
used the density PS fro m other redsh i ft sur veys in es- 
timating (3. Recently, Zaroubi et al. (2002) have ap- 
plied UMV estimator to SEcat, a combination of the 
SFI and the ENEAR catalogs, and found consistent /3 
estimates from both S-S and v-v comparison methods. 
Interestin gly, low (3 est imation (/3 = 0.2-0.4) has b een re- 
ported bv lRiess et aT] l|1997|) and IBlakeslee et alJ (1999) 
who have compared the Optical Redshift Survey (ORS; 
iSantiaeo et al.lll9'95jl with peculiar velocities from Type 
la supernovae and a surface brightness fluctuation (SBF) 
survey of galaxy distances, respectively. 

Our measurement of the mass fluctuation level 
er^m 6 = 0-56io'2i is larger than that derived from 
the Wilkinson Mic rowave Anisotr opy Probe (WMAP ; 
IBennett et aTll2003j1 data analysis. ISpergel et all ( 2003) 
have used WMAP data to derive a value of cgf^ 6 = 
0.44 ± 0.10. The reason for the large cs^™ 6 is that our 
measurement of the rms fluctuation of galaxy distribu- 
tion gives somewhat large and uncertain value cis.s = 
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TABLE 4 

Comparison of j3 and us^m 6 Derived from the Other Peculiar Velocity Samples 



Parameter 


Data 




Method Ref. 


P 


0.50 ±0.10 


ENEAR, IRAS PSCz 




OME (v-v) 1 


u ' oi -0.06 


SEcat, IRAS PSCz 




UMV (v-v) 2 


n 57+° 11 
U -°'-0.13 


SEcat, IRAS PSCz 




UMV (S-S) 2 


0.50 ± 0.04 


Mark III, IRAS 1.2 Jy 




VELMOD (v-v) 3 


1.20 ±0.10 


Mark III, IRAS 1.2 Jy 




POTENT (S-S) 4 


0.89 ±0.12 


Mark III, IRAS 1.2 Jy 




POTENT (5-8) 5 


0.50 ±0.06 


D n -u, IRTF, ESO, UGC 




X 2 fit. (v-v) 6 


0.39 ±0.17 


SMAC, IRAS PSCz 




X 2 fit. (v-v) 7 


0.40 ±0.15 


SNIa, IRAS 1.2Jy 




X 2 fit. (v-v) 8 


0.30 ±0.10 


SNIa, ORS 




X 2 fit. (v-v) 8 


0.55 ±0.06 


SNIa, IRAS PSCz 




X 2 fit. (v-v) 9 


42+ - 10 


SBF, IRAS 1.2 Jy 




X 2 fit. (v-v) 10 


0.26 ±0.08 


SBF, ORS 




X 2 fit. (v-v) 10 


a 8 n° m 6 


50+ ' 25 
U - DU -0.14 


SCI 




Vrms 1 1 


0.5ll° o| for T = 0.25 


ENEAR 




CF 12 


1 ' ± -0.35 


ENEAR 




ML 13 


0.88 ±0.15 


Mark III 




ML 14 


0.49 ±0.06 


Mark III 




ML 15 


(0.71-0.77)±0.12 


Mark III 




POTENT 4 


References. (1) INusser et all (1200111: ( 21 IZaroubi et al J (1200211 : f 31 
IWillick & Straus^ (l998|k (4l|Kolatt & Dekel| l|l997fl : (5l|Sisad et alj(1998f: (6l|Hudsorj 
1 199411: f7l~lHudson et alt 12004); (8) IRiess et alJ (1199711: C91 IRadburn-Smith et al. 
120011: (101 IBlakeslee et alJ (|1999Tl: fill IBoreani et all \199W (12) IBoreani et al 


l^QOOT): C131 IZaroubi et alJ 12001D: C141 IZaroubi et alJ 11997D: (151 ISilberman et al. 





1.14^0 34- By assuming that WMAP cosmological pa- 
rameters cr s = 0.9 and f2 m — 0.29, we obtain cr^g = 

asfl^/Ps = 0-87to.ll for #3 = 0.489. Despite the 
higher amplitude of mass fluctuation than the WMAP 
value, our estimation of the matter density PS Sl^P$(k) 
is lower than those of other studies (Table 1 and Fig. 9). 
For example, at k = 0.1047 /iMpc" 1 (A = 60 h^Mpc), 
our PS value is Vl^P s {k) = 2.blt° f 4 x 10 3 (h^Mpc) 3 , 
which is lower compared with the high-amplitude power 
spectra found from the ML analyses of all-sky peculiar 
velocity samples such as Mark III, SFI, and ENEAR. 
At k — 0.1 /iMpc -1 scale, the matter density powers 
(0,^ 2 P S ) f rom the ML analys es are (4.8 ± 1.5) x 10 3 
/Mark III: IZaroubi et a!Jll997ll . (4.4 ± 1.7) x 10 3 (SFI; 
iFreudline et alJll999D. "and (6 .5 ± 3) x 10 3 (/i^Mpc) 3 
(ENEAR; Zar oubi et al.H200l|l . Those analyses also give 
the high mass fluctuation levels (crs^™ 6 ) around .8-1.1 
(Table 3 and 4) . On the contrarv. lSilberman et all (|200lH 
found a power deficiency at k = 0.1 ftMpc -1 from the 
ML analysis of the SFI sample, and their density PS and 
(Js^nf are consistent with our results. 

As shown in Figure 8, the f3 in the quasi-linear regime 
is more accurate than that in the linear-regime. This 
implies that if we rely on the linear regime where the 



density and velocity PS are simply related by the lin- 
ear gravitational instability theory as in equation (5), 
we need peculiar velocity data with large survey volume 
for accurate determination of the (3 parameter. On the 
other hand, if we use the full information of the mo- 
mentum field from linear to quasi-linear regimes where 
the relation between the density and momentum PS are 
accurately known, a peculiar velocity sample gives more 
accurate (3 compared with that estimated only in the lin- 
ear regime. For the best (3 determination, it is essential 
to derive a more accurate relation that connects the den- 
sity PS with the momentum PS beyond the limit of the 
quasi-linear regime by applying the higher order pertur- 
bation theory. 
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